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Using the (2+l)-dimensional viscous hydrodynamic code VISH2+1 with a temperature dependent 
specific shear viscosity (ry/s)(T), we present a detailed study of the influence of a large hadronic shear 
viscosity and its corresponding relaxation time r,r on the transverse momentum spectra and elliptic 
flow of hadrons produced in 200 A GeV Au+Au collisions. Although theory, in principle, predicts 
a well-defined relation t-^T — it{T) x (rj/s) (T), the precise form of k(T) for the matter created in 
relativistic heavy-ion collisions is not known. For the popular choice k = 3 the hadron spectra are 
found to be insensitive to a significant rise of r^/s in the hadronic stage, whereas their differential 
elliptic flow V2{pt) is strongly suppressed by large hadronic viscosity. The large viscous effects 
on V2 are strongly reduced if (as theoretically expected) k{T) is allowed to grow with decreasing 
temperature in the hadronic stage. This implies that, until reliable calculations of ^(r) become 
available, an extraction of the hadronic shear viscosity from a comparison between VISH2+1 and a 
microscopic hadron cascade or experimental data requires a simultaneous fit of {ri/s){T) and k(T). 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld, 24.10.Nz 
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I. INTRODUCTION 

A fluid state of matter, quark-gluon plasma (QGP), 
is created in ultra-relativistic heavy-ion collision experi- 
ments at the Relativistic Heavy Ion Collider (RHIC) [1- 
4] . Theoretical analysis of these experiments established 
that QGP behaves like an almost perfect liquid with very 
small viscosity [5-7] . Much effort has been focused on de- 
termining the QGP transport parameters, in particular 
its specific shear viscosity (77/s)qgp, he. the ratio be- 
tween its shear viscosity ry and entropy density s (see 
[8, 9] for recent reviews). 

In a recent article [10] the newly developed hybrid 
code VISHNU [11] has been used to extract (??/s)qgp 
from the observed collision centrality dependence of the 
integrated charged hadron elliptic flow V2. This code 
couples the macroscopic evolution of the QGP by the 
(2-|-l)-dimensional viscous hydrodynamic code VISH2+1 
to the Boltzmann cascade UrQMD which describes the fi- 
nal hadronic rescattering and freeze-out stage microscop- 
ically. The microscopic simulation of the late hadronic 
stage is numerically costly, and a macroscopic descrip- 
tion with viscous fluid dynamics would therefore be much 
preferred if valid. Unfortunately, a detailed study pre- 
sented in [11] indicated that the microscopic UrQMD dy- 
namics cannot be faithfully simulated with viscous hy- 
drodynamics if one assumes the frequently used relation- 
ship T^r = 3-j between the specific shear viscosity r]/s 
and the microscopic relaxation time r^r for the shear vis- 
cous pressure tensor tt'"^, scaled with the temperature 
T of the fluid. Relations of the type Tt^T=k,^, with 
constant proportionality factors k, are found theoreti- 
cally in both the extreme weak-coupling (for a masslcss 
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Boltzmann gas one finds k = 6 in Israel- Stewart theory 
[12, 13] and k = 5 in the modified approach by Denicol 
et al. [14]) and extreme strong-coupling limits (where 
one has k = 4-2 In 2 = 2.614 for A/" = 4 super- Yang-Mills 
theory at infinite coupling [15-17]). Other recent work, 
however, based on modified Kubo relations [18-21] and 
a deeper analysis of the Boltzmann equation and its con- 
nection to viscous hydrodynamics [14, 22, 23], suggests 
(in some cases strong) temperature dependence of k. 

Here we will explore one such proposed relation, 
K = {e+p)/p [19] (where e and p are the energy density 
and pressure of the system), which leads to a strong in- 
crease of k(T) with decreasing temperature in the mas- 
sive hadron resonance gas below the quark confinement 
temperature Tc- Such an increase is qualitatively consis- 
tent with certain observations made in the recent VISHNU 
study [11]. We here use VISH2+1 to investigate, within a 
purely hydrodynamic framework, systematically the con- 
sequences of increasing shear viscosity and shear pres- 
sure relaxation time in the late hadronic stage on the 
transverse momentum spectra and elliptic fiow of soft 
(pt < 2 GcV/c) hadrons produced in Au+Au collisions at 
RHIC. Our work differs from an earlier study by Bozek 
[24] of the effects of temperature-dependent specific bulk 
and shear viscosities by focusing on shear viscosity and 
investigating situations in which the shear viscosity of the 
hadron gas is larger than that of the QGP (rather than 
the other way around [24]), as expected on basic theo- 
retical grounds [25]. While the present work was being 
completed, a related study appeared [26] which focusscs 
chiefly on the question whether recent data from Pb-|-Pb 
collisions at the Large Hadron Collider (LHC) [27] re- 
quire an increase of the QGP shear viscosity with rising 
temperature. 

The paper is organized as follows: In Sec. II we briefly 
review the viscous hydrodynamic model and discuss the 
specific ingredients used in the present study. The ef- 
fects of a large hadronic specific shear viscosity (77/s)hg 
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on the fireball evolution are discussed in Sec. III. In 
Sec. IV we discuss the dependence of the transverse mo- 
mentum spectra and elliptic flow of emitted hadrons in 
Au+Au collisions on (77/s)hG; the decoupling tempera- 
ture Tdoc, and the collision centrality. Section V is dedi- 
cated to a detailed discussion of the viscous corrections to 
the freeze-out phase-space distribution and their effects 
on spectra and elliptic flow. All results up to this point 
assume a constant factor k = 3 in the relation t^^T ~ 
between the specific shear viscosity 77/s and the micro- 
scopic relaxation time Tt^; in Sec. VI we explore the con- 
sequences of making k{T) temperature dependent and 
letting it grow during the quark-hadron phase transition. 
A final discussion in Sec. VII concludes our paper. 



II. VISCOUS HYDRODYNAMICS: SPECIFIC 
INGREDIENTS FOR THE PRESENT STUDY 

VISH2+1 [28] solves the second-order Israel-Stewart 
equations for causal relativistic viscous fluid dynamics 
[12] in the spatial plane transverse to the beam direction 
and in time, assuming boost-invariance of the longitudi- 
nal expansion. To avoid repetition we refer the reader 
interested in the technical details to earlier descriptions 
of the specific form of the evolution equations and the 
equation of state s95p-PCE used here (see specifically 
Sections II and III in Ref. [29]). The energy- momentum 
tensor is decomposed as T^"^ = eu'^u'^ —pA^'^ +tt'^'^ where 
tt'^'^ is the viscous pressure tensor. We consider only 
shear viscosity and ignore bulk viscous effects; in this 
situation the Israel-Stewart equations describe the evo- 
lution of TT^'^ towards it Navier-Stokes limit 2rja^^ on 
a microscopic relaxation time scale T^r, where 77 is the 
shear viscosity and a^'^ is the velocity shear tensor which 
evolves hydrodynamically in space and time. We initial- 
ize TT^"^ with its Navier-Stokes value tt^"^ — '^V'^q'^ at ini- 
tial time To, calculated from the initial velocity profile 
u'' = (m^, u^, uy, u'') = (1,0, 0, 0). 

The generation of hydrodynamic fiow from the pres- 
sure gradients in the system is controlled by the 
fluid's equation of state (EOS) for which we use s95p- 
PCE [29, 30] with chemical decoupling temperature 
Tchcm = 165 MeV. This EOS interpolates between state- 
of-the-art Lattice QCD data at high temperatures and a 
chemically frozen hadron resonance gas at low temper- 
atures. Chemical freeze-out at Tchem = 165 MeV guar- 
antees that the final hadron yields, calculated by in- 
tegrating the final hadron momentum spectra obtained 
from the hydrodynamic output along an isothermal de- 
coupling surface of temperature Tdcc < Tchcm via the 
Cooper-Frye procedure [31] followed by resonance de- 
cay [32, 33], agree with experimental measurements in 
200AGeV Au-^Au cohisions at RHIC [3, 34, 35]. At de- 
coupling, we parametrize the local distribution function 
in the Cooper-Frye formula by a local thermal equilib- 
rium function plus a small viscous correction which de- 
pends on the value of the viscous pressure tensor tt'^'^ on 



the freeze-out surface and increases quadratically with 
particle momentum [29, 36]. Unless noted otherwise, we 
use Tdcc = 120 MeV. 

We initialize the hydrodynamic evolution with an en- 
ergy density profile obtained from the optical f KLN model 
[37-39]. The model yields the initial gluon density dis- 
tribution which, after thermalization, gives directly the 
initial entropy density which is then converted to energy 
density using the EOS s95p-PCE. The normalization of 
the initial entropy density is adjusted in the most cen- 
tral collisions to reproduce the finally measured charged 
hadron multiplicity. Due to viscous entropy production, 
changing 77/s requires a readjustment of this normaliza- 
tion to keep the final multiplicity fixed. After normal- 
ization in central collisions, the centrality dependence of 
the final charged hadron multiplicity is obtained directly 
from the f KLN model, without further adjustment of pa- 
rameters. 

The key ingredients whose influence on the genera- 
tion of radial and elliptic flow we want to study here 
are the temperature dependence of the specific shear vis- 
cosity r]/s and of the proportionality constant between 
rj/s and the temperature-scaled microscopic relaxation 
time TttT, K=^^yj. Specifically, we will explore scenar- 
ios where r]/ s = 0.16 is a constant in the QGP phase but 
increases by variable amounts during the transition from 
QGP to hadrons, using the following parametrization for 
its temperature dependence: 



■>l^rj.^ ^ (??/g)QGP + (?7/g)HG 

+ (V^)qgp - (V^)hg ^^^^^ ( ^J-T. 



Tr 
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Here Tc = 170MeV, and (77/s)qgp = 0.16 and {r]/s)B.G 
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FIG. 1: (Color online) Five choices for the temperature de- 
pendent {rj/s){T) studied in this work. The constant values 
at low T are multiples of 0.08 ~ ^i^. 

arc (different) constants for the QGP and HG (hadron 
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gas) phases. We will explore the range 0.16 < (ry/s)HG 
<0.48, as illustrated in Figure 1. In the next three sec- 
tions K will be held constant at k = 3;^ consequences of a 
temperature dependent k[T) = {e+p)/p will be explored 
in Sec. VI. 



III. HYDRODYNAMIC EVOLUTION 

In order to study how the fireball evolves with a tem- 
perature dependent {r]/s){T) that increases in the HG 
phase, we graph the time evolution for the average trans- 
verse flow velocity {{vj_)) (the average over the transverse 
plane being defined with the lab-frame energy density 



7^e as weight), the spatial eccentricity e^, — ^^^2 ^^2^^ 
the lab-frame energy density distribution, the flow mo- 



ot 



mentum anisotropy Ep 



(where (...) denotes 



simple integration over the transverse plane and Tq'^ is 
the ideal fluid part of the energy- momentum tensor, with- 
out viscous pressure contributions) , and the total momen- 
turn anisotropy e'p = ^y^^^^y^^ for different choices of the 
temperature dependence of ry/s. 

Since shear viscosity leads to viscous heating which 
generates entropy, holding the finally observed hadron 
multiphcity fixed requires that an increase in {r]/s){T) 
must be accompanied by a decrease of the initial entropy 
of the fireball. We implement this by a decrease of the 
normalization of the initial entropy density distribution, 
keeping its shape fixed. Whereas for fixed initial condi- 
tions an overall increase oirj/s leads to stronger radial ac- 
celeration due to a positive contribution from the viscous 
pressure tensor tt'"^ to the transverse pressure gradients 
[36, 40-43], this effect is largely compensated [26, 29, 44] 
after rescaling the initial entropy density to ensure fixed 
final multiplicity. For our temperature-dependent rj/s 
this compensation no longer works in the same way: af- 
ter rescaling the initial entropy density profile, to com- 
pensate for increased viscous heating in the hadronic 
phase, the QGP core shrinks and the HG corona grows 
in size. Since the viscous pressure is relatively larger in 
the hadronic phase than in the QGP, the effective trans- 
verse pressure gradient is reduced when increasing rj/s 
only in the hadronic phase, leading to weaker radial ac- 
celeration. This can be seen in Fig. la, where we see a 
reduction oof the growth rate of the average radial flow 
velocity {{v±)) with increasing values of (?7/s)hGi holding 
(??/s)qgp = 0.16 fixed. 

The larger shear viscosity in the hadronic corona leads 
to a more rapid initial decay^ of the spatial fireball eccen- 
tricity Sx (see Fig. 2b) and a slower growth rate and lower 
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FIG. 2: (Color online) The average radial flow {{i'±}), spatial 
eccentricity Sx , and the flow and total momentum anisotropies 
Ep and e'p for Au+ Au collisions at 6 = 7.5 fm as functions of hy- 
drodynamic evolution time r— ro, for tq = 0.4 fm/c and kinetic 
freeze-out temperature Tdcc^lSOMeV. Lines with different 
symbols correspond to different temperature dependences of 
rj/s as shown in Fig. 1. 



asymptotic value of the flow momentum anisotropy Ep 
(Fig. 2c, open symbols). The spatial eccentricity curves 
in Fig. 2b all cross around t—tq = 4.5 fm/c, indicating 
the transition from stronger decay of Sx at early times to 
weaker decay at late times for larger values of (?7/s)hg- 
This is a consequence of the reduced flow anisotropy Ep 



^ The specific values (»;/s)qgp =0.16 and k = 3 diosen here agree 
with those used by us in the earlier studies [10, 29] while the 
recent work [26] assumes re = 5. 

^ We note that e is defined by integrating at fixed time r over the 
entire transverse plane, including both thermalized and already 



decoupled matter. It is possible that the strong initial decay of 
e seen in Fig. 2b arises mostly from contributions in that part of 
the hadronic corona that has already decoupled. 
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FIG. 3: (Color online) (a): Transverse momentum spectra for charged hadrons, pions, and protons from VISH2+1 for the 5% 
most central Au+Au collisions (6 = 2.33fm). (b-c): Differential elliptic flow V2{pt) for charged hadrons (b), pions (c) and 
protons (d) from Au+Au collisions at 20-30% centrality (6 = 7.49fm). Lines with different symbols correspond to different 
values of (77/s)hg as shown in Fig. f; Tdcc = 120 MeV. Decay products from all strong resonance decays are included. Charged 



hadrons include tt , K , p, S 



f2 , and their antiparticles. 



shown in Fig. 2c. 

The lines with filled symbols in Figure 2c show that 
the effects of increased hadronic viscosity on the asymp- 
totic values of the total momentum anisotropy are 
much stronger than on the flow anisotropy Sp: while 
the latter decreases by about 25% from (ry/s)HG = 0.16 
to (7?/s)hg =0.48, the corresponding decrease for e'p is 
almost twice as large. Also, most of the effect on e'p 
happens at late times t—tq > 4.5fm/c when most of the 
matter has converted into hadron gas. This reflects the 
growth of the Navier-Stokes value Trfjg = 2r](j^^'^ of the vis- 
cous pressure contribution to T^'^ in the hadronic phase 



where Tync increases. In contrast to Sp, the total momen- 
tum anisotropy e'p does not saturate at late times after 
the spatial eccentricity (which drives the flow anisotropy) 
has essentially decayed to zero; its continued increase is 
due to the continuing decrease of the magnitude of the 
tt'^" components whose contribution to e'p is negative [43] . 

The large difference between the late-time values of 
ep and e'p for high values of 77HG shows that, for strong 
hadronic viscosity, the viscous corrections to the local 
thermal equilibrium distribution on the kinetic decou- 
pling surface at Tdec are big. We will explore this in 
more detail in Sec. V. 
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FIG. 4: (Color online) Similar to Fig. 3, for fixed (?7/s)hg =0.48 and different decoupling temperatures Tdcc ranging from 160 
to 100 MeV. In panels (b)-(d), dotted lines show V2{pt) calculated without the Sf correction whereas the solid lines show the 
full calculations. 



IV. SPECTRA AND ELLIPTIC FLOW 



A. Central and semi-peripheral Au+Au collisions 



Figure 3 shows the transverse momentum spectra for 
charged hadrons, pions and protons from central Au+Au 
collisions (0—5% centrality) and their elliptic flows V2{pt) 
for semiperipheral Au+Au collisions (20—30% central- 
ity) for different choices of the hadronic shear viscosity 
(??/s)hg- The pT-spectra in Fig. 3a are seen to be com- 
pletely insensitive to the value of (77/s)hg- From the 
reduction of the radial flow seen in Fig. 2a one would 
have expected steeper spectra for larger (77/s)hg since 
Tdoc = 120MeV is held fixed; clearly, for pr < 2GeV/c, 



the viscous correction 6f to the local equilibrium distri- 
bution at freeze-out (which will be analyzed in greater 
depth in Sec. V) happens to almost exactly compensate 
for the loss of radial flow, over the entire range of (?;/s)hg 
values studied here. 

This is not true for the elliptic flow which is 
strongly reduced when the hadronic viscosity is increased 
(Figs. 3b-d). For protons a striking effect is seen for 
('7/'S)hg > 0.32: The proton elliptic flow turns negative 
(i.e. protons show stronger flow perpendicular than par- 
allel to the reaction plane) for low px- This effect is 
caused entirely by the Sf correction. Sf grows not only 
with pt, as is well known, but also with the mass of the 
hadron. For massive hadrons, the shear viscous Sf cor- 
rection can be a strong effect even at = 0. In Fig. 3 
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negative V2(j>t) caused by Sf at low px is not visible for 
pions, but for protons and would be much stronger for 
hyperons or J/ tp mesons [45] if they also followed viscous 
hydrodynamical evolution down to T^ec = 120 MeV. 

The effect of the 6f correction is studied in Fig. 4, for 
various choices of the decoupling temperature T^oc ■ We 
hold the hadronic shear viscosity fixed at (?7/s)hg = 0.48, 
the largest value studied here. The effect of variations 
in Tdoc on the spectra in Fig. 4a is similar to what we 
observed in [29]: lower decoupling temperatures cause 
flatter proton spectra due to larger radial flow, steeper 
pion spectra due to the cooling effect which dominates 
for light particles, and almost no change in the charged 
hadron spectra whose mix of light and heavy particles 
effectively balances the counteracting cooling and radial 
flow effects. 

In Figs. 4b-d we plot the differential elliptic flow for 
charged hadrons, pions and protons. The dotted lines 
show a calculation that ignores the viscous Sf correction 
at freeze-out and thus only includes the Tdec-dependence 
of the pure flow effects. We see that lower Tdcc values 
suppress V2{pt) for protons but increase it for pions at 
low pt- This is really a consequence of the accompanying 
change of the p^-spectra: Due to the large hadronic vis- 
cosity, very little additional flow momentum anisotropy 
is generated at temperatures below Tc. However, due 
to cooling, the pion spectra get steeper with decreas- 
ing Tdoc, moving more of their momentum anisotropy 
to low transverse momenta which leads to the increase 
of pion V2{pt) at low pt- Conversely, the proton spectra 
get flatter, in spite of cooling, due to additional radial 
flow developing between Tc and Tdcc! consequently, their 
total momentum anisotropy gets shifted on average to 
larger transverse momenta, causing a reduction of pro- 
ton V2{pt) at low Pt (accompanied by an increase at 
high Pt < 2 GeV/c, beyond the range shown here). Both 
the flattening of the proton spectra and the shifting of 
their elliptic flow to larger pT are stronger for the case 
of large hadronic shear viscosity ((ry/s)HG = 0.48) stud- 
ied here than for the case of temperature-independent 
r]/s = 0.16 studied in [29]: The large hadronic viscosity 
generates stronger additional radial flow but less addi- 
tional momentum anisotropy in the hadronic stage than 
does constant 77/3 = 0.16. Note that, without 5f, proton 
V2{pt) never turns negative, even for the largest hadronic 
shear viscosity studied in this work. 

The solid lines in Figs. 4b-d show the full calculation of 
V2{Pt) including the 5f correction. We see larger Sf ef- 
fects for protons than pions, due to their larger rest mass 
[43] . The full calculations feature a non-monotonic vari- 
ation of pion and charged hadron V2{pt) with decoupling 
temperature Tdcc: The suppression from Sf is smaller for 
Tdoc = 160 MeV than for Tdcc 140 MeV. The like expla- 
nation is that Tdec = 160 MeV is so close to the inflec- 
tion point Tc of the shear viscosity {rj/ s){T) that, due to 
the finite relaxation time r^r ~ 2fm/c at this tempera- 
ture, the viscous pressure tensor has not yet had time to 
fully evolve to its (larger) hadronic Navier-Stokes value 



whereas at Tdoc complete relaxation has been achieved. 
At sufficiently low Tdec, Sf decreases with decreasing the 
decoupling temperature, since now rj/s has reached its 
new, higher hadronic level and -k^" becomes smaller sim- 
ply due to hydrodynamic expansion [43]. 



B. Minimum bias collisions 

In Figure 5 we show pT spectra and differential elliptic 
fiow for charged hadrons, pions and protons from min- 




P, (GeV) 

FIG. 5: (Color online) Transverse momentum spectra (a) and 
differential elliptic flow V2{pt) for charged hadrons (b) and 
pions and protons (c) from minimum bias Au+Au collisions 
at RHIC, for various {r]/s)(T) as indicated (c. /. Fig. 1). 
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imum bias Au+Au collisions with T^ec = 120 MeV. For 
these we summed our calculated results over all collision 
centralities < 80%. The dependence on collision central- 
ity is discussed in the next subsection. 

Similar to what we saw in Fig. 3a for central Au+Au 
collisions, the spectra shown in Fig. 5a exhibit almost no 
sensitivity at all to variations of the specific shear vis- 
cosity (77/s)hg hi the hadron gas stage. We did observe 
some flattening of the charged hadron spectrum in the 
most peripheral (70—80%) centrality bin studied, where 
the viscous effects are strongest and the Sf correction is 
largest. Due to its low weight in the average, this weak 
effect is not visible in the minimum bias result. 

In Figures 5b and c, the minimum bias differential 
V2{pt) of all charged hadrons, pions and protons are 
shown for different (77/s)hg- We see that the features 
observed in Fig. 3 for the specific 20—30% centrality 
bin carry over, qualitatively unchanged, to event sam- 
ples without centrality selection: a significant increase 
of rj/s in the hadron gas phase has a strong suppres- 
sion effect on V2{pt)- However, as shown in Sec. IV A, 
the suppression arises mostly from the 5f correction at 
kinetic freeze-out, with a much smaller contribution ac- 
counting for the lack of growth of the total momentum 
anisotropy in the hadronic phase when (?7/s)hg becomes 
large. Hence, the strong suppression of differential ellip- 
tic flow by large hadronic shear viscosity shown here is 
critically dependent on the validity of viscous hydrody- 
namics as the correct framework for evolving 5f all the 
way down to Tdcc = 120 McV. This is assumed here, but 
not supported by the analysis presented in [II]. 
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FIG. 6: (Color online) Eccentricity-scaled charged hadron 
elliptic flow V2/£ as a function of the multiplicity density 
{l/S){dNch/dy), for different values of {rj/s)HG- The over- 
lap area S — tta/ {{x^}) {{y^}) is calculated from the same initial 
profiles as the spatial eccentricity e. 



V2{pt) around pr = 0.5 GeV/c (i.e. close to their average 
Pt), caused by large 5f corrections at freeze-out.'^ 



C. Centrality dependence of elliptic flow 



V. Sf CONTRIBUTIONS 



The centrality dependence of the eccentricity-scaled el- 
liptic flow V2/e is shown in Fig. 6 where we graph this 
quantity as a function of the final charged multiplicity 
density {1/ S)dNch/dy. for different values of (77/s)hg. 
(We obtain V2 by integrating V2{pt) over all pr, with- 
out regard to possible pt cuts imposed by experimental 
constraints.) Strong suppression of V2/ s by hadronic vis- 
cosity is observed even in the most central collisions, but 
the effect is stronger in peripheral collisions. An increase 
of (77/s)hg thus not only decreases i'2/e, but also changes 
the slope of its centrality dependence. We note in pass- 
ing that in recent studies with the hydro-|-cascade hybrid 
code VISHNU [10] this slope was fixed and controlled by 
the effective dissipation encoded in the hadron cascade, 
and that in [11] an (unsuccessful) attempt was made to 
extract the temperature-dependence of (?//s)hg (here as- 
sumed to be T-indcpendent) by matching the magnitude 
and slope of the corresponding V2/e vs. {l/S)dNch/dy 
curves from VISH2+1 to those from VISHNU. We also ob- 
serve that for the largest value of (77/s)hg studied here, 
(??/s)hg =0.48, the total charged hadron elliptic flow 
turns negative in the most peripheral (70—80%) central- 
ity bin. We found that this is caused by negative pion 



Due to non-zero viscous pressure components tt^'^ , the 
distribution function fi(x,p) for hadron species i must 
deviate on the freeze-out surface from local equilibrium: 



fl{x,p) = fcq,i{x,p) + Sfi{x,p). 



(2) 



We use [36] 



Sfi = /c. 



2 T2 e + p' 



(3) 



3 For 70-80% centrality and (r?/s)HG = 0.32-0.48, we found for 
that for pions V2{pt) first rises at very low pT , then turns nega- 
tive for 0.25 <Pt < 0.75 GeV/c before turning positive again and 
continuing to grow approximately linearly with px- This is dif- 
ferent from protons whose V2{pt) turns negative right away at 
small px, again with a minimum around 0.5GeV/c. All these 
effects are caused by large Sf effects; in this centrality bin we do 
not trust viscous hydrodynamic predictions to be very robust. 
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(GeV) (GeV) p^ (GeV) 

FIG. 7: (Color online) The Sf correction for constant 77/5 = 0.16 to the pr-spectra for charged hadrons (a), pions (b) and 
protons (c) at 0—5% centrality, and to their differential elliptic flow (d-f) at 20—30% centrality. Lines with different symbols 
denote individual contributions as described in the text. 



noting that also other forms have been suggested in the 
hterature [46, 47]. The numerator can be written as 



= 7r''^(x) 
+ S(a;) 
+ A(x) 



^|(2cosh^(,-,)-l) + 4^ 

sm(20„) 

oPT sin(0p+(/)^,) 
-2 — rriT cosh(y-7?) — . ^ 

2 • , 2. \ , Pt f. sin(20p) 
-mj. smh (y-77) + — I 1 - . 

2 V sm(2(p„) 



sin((/)p-0^) 
tan(20„) 



PxiTiT cosh(y— 77)7 

p"^ sin{2{4)p~4>y)) 
2 sin(2(/.„) 



sm 



sin(20^,) 



(4) 



where S = tt^^+tt?'?', A^tt-^^- 



Because of boost- 



invariance, tracelessness and orthogonality to it'^, only 
three components of tt^'^ are independent; we take them 



as E, A, and tt"^"^. mT — yrn? 



■p^ is the transverse 



mass of the particles, 0p is the azimuthal angle of pr, 
and 4>v{x) is the azimuthal angle of the fluid velocity v 
at point X. 

We now discuss the individual contributions from 
Eq. (4) to the pT-spectra and elliptic flow, for the 
cases of constant 77/s = 0.16 (Figure 7) and temperature- 
dependent {t]/s){T) (Figure 8). In panels (a-c) we 
show the fractional contribution SN/Ncq from Sf to the 
Cooper-Frye spectra of charged hadrons (a), pions (b) 
and protons (c). At low px, the contributions propor- 
tional to tt'^'^ and A (first and last terms on the r.h.s. 
of Eq. (4)) are small and overshadowed by the contri- 
bution from the average transverse viscous pressure E. 
The first (negative) term ~ —m^ in the expression mul- 




FIG. 8: (Color online) Similar to Fig. 7, but for temperature dependent ('7/s)(T), Eq. (1), with (?7/s)hg =0.48. 



tiplying S dominates at low pT ■ It obviously grows with 
rest mass, leading to large negative 6N/Neq corrections 
at low pt for heavy hadrons such as and J/ ip [45] . For 
protons the effect remains below 10% in central Au+Au 
collisions, i.e. Sf corrections arc small and the calcula- 
tion is reliable. At larger px, all three contributions in 
Eq. (4) turn positive and 5N/Ncq switches sign (around 
0.5GcV/c for pions and around IGcV/c for protons). 
Again, the term ^ E first dominates, but since it grows 
only linearly at large pr it is eventually (at pr^^ GeV/c) 
overtaken by the term ^tt"^"^. For constant Ty/s^O.lB, 
\dN/Ncq\ remains below 25% up to px = 2GeV/c for all 
three spectra shown, ^ and the calculation is therefore reli- 
able. For large hadronic viscosity (?7/s)hg = 0-48 (Fig. 8) 



* The Sf effects on charged hadron spectra can be qualitatively 
understood from those on pion and proton spectra by noting that 
at low Px charged hadrons are dominated by pions whereas at 
larger px heavier hadrons become increasingly more important. 



the Sf corrections to the pr-spectra are larger, in partic- 
ular the term ~7r'''", and \SN/Ncq\ reaches 70-80% at 
Pt = 2 GeV/c, indicating the imminent breakdown of the 
viscous hydrodynamic expansion ^ /cq. 

In the lower panels of Figs. 7 and 8 we show the 
6f contributions to the differential V2(j>x) for charged 
hadrons (d), pions (e), and protons (f), again separated 
into their individual contributions according to Eq. (3). 
We see that for low px all three terms in Eq. (4) con- 
tribute to the suppression of elliptic flow, but that in 
this case at high px the term proportional to the vis- 
cous pressure anisotropy A = 7r^^— 77**^ plays the domi- 
nant role, overshadowing the terms S and (except for 
the largest hadronic viscosities) also ^tt'^'^. The latter 
grows quadratically with px and eventually wins over the 
term S; for large hadronic viscosity (Fig. 8) it even ex- 
ceeds the anisotropy term ~ A at sufficiently large px- 
The term proportional to the average transverse viscous 
pressure S individually generates a positive elliptic flow 
correction at large px (i.e. at pT ^ 2 GeV/c for constant 
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r^/s = 0.16 and at ^ (1 — 1-5) GcV/c for T-dependcnt 
{r]/s){T) with (77/s)hg = 0-48). Similarly the anisotropy 
term ^ A by itself increases proton elliptic flow at low 
Pt if the hadronic viscosity is large enough (Fig. 8f). In 
the sum, however, these positive individual corrections 
are always overwhelmed by the remaining two negative 
corrections, leading to an overall suppression of V2{pt) 
at all Pt in all cases. Interestingly, the negative proton 
elliptic flow at low pt and large (r?/s)HG values noted 
earlier (Figs. 3-5) is not caused by the viscous pressure 
anisotropy A, but by the average transverse viscous pres- 
sure E (green triangles in Fig. 8f). This phenomenon is 
driven by the eS'ect of E on the proton spectra (Figs. 7c 
and 8c): S suppresses the spectra at low pT^ leading (in 
extreme situations) to the formation of a shoulder in the 
proton spectra which is known [48] to cause negative V2- 

VI. LARGE HADRONIC RELAXATION TIMES 



viscosity is kept at (?7/s)qgp = 0.16 throughout, but we 
toggle (?7/s)hg in Eq- (1) between the two values 0.16 
and 0.48 (see Fig. 1). 

Figure 10 shows a similar analysis as Fig. 2, but now 
comparing constant with T-dependent k values. From 
Fig. 10a wc conclude that the temperature dependence 
of K has no visible influence on the evolution of the spa- 
tial eccentricity Sx, irrespective of whether the specific 
shear viscosity rj/s grows in the hadronic phase or not. 
On the other hand we see in Fig. 10b that a n{T) that 
grows around and below as shown in Fig. 9 reduces sig- 
nificantly the viscous suppression of the total momentum 
anisotropy that is otherwise caused by a large hadronic 
shear viscosity.^ Analyzing panel (b) of Fig. 10 in more 
detail, we observe that during the early stage of the evo- 
lution larger hadronic relaxation times have little effect 
on the flow momentum anisotropy £p, consistent with the 
almost unchanged decay rate of the spatial eccentricity 



Motivated by the study of the VISHNU model in [11] 
we explore in this section the consequences of very large 
relaxation times in the hadronic phase. Specifically, 
we assume a relation proposed in [19], 



n{T) 



e+ p 
P 



(5) 



which can be easily worked out for our EOS s95p-PCE 
and is shown in Fig. 9. In the massless limit (i.e. at large 
T where the EOS approaches e = 3p), this expression ap- 
proaches the value k = 4. To explore effects specifically 
related to the T-dependence of k, we compare in this 
section results from Eq. (5) with those for constant k = 4 
(and not k — 3 as in the preceding sections). The QGP 




0.2 0.3 0.4 
T (GeV) 



0.6 



FIG. 9: The temperature dependent n{T) from Eq. (5) for 
EOS s95p-PCE (solid), compared with the massless limit 
K = 4 (dashed). 
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FIG. 10: (Color online) Same as Fig. 2b,c, but for k{T) from 
Eq. (5) and constant k = 4 instead of ac = 3. For the QGP 
(77/s)qgp =0.16 is used throughout whereas {ri/s)iiG is varied 
between 0.16 and 0.48 as indicated in the legend. 



Please note that the extremely rapid rise of k(T) below T ~ 
50 MeV seen in Fig. 9 is irrelevant in this context because the 
fireball matter decouples already at T^^^ = 120 MeV. 
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seen in panel (a) that drives the anisotropic flow. At late 
times, however, the larger k{T) is seen to have a small 
positive effect on the generation of anisotropic collective 
flow. Increasing the response time t-^ with which the 
viscous pressure tensor tt'^" can react to changes in the 
velocity shear tensor apparantly allows the collective flow 
anisotropy to grow more easily, with less viscous damp- 
ing, than if tt^'' is allowed to relax to its Navier-Stokes 
value Trfjg = 2ria^'^ more quickly. This is a cumulative ef- 
fect that becomes visible most clearly at late times when 
most of the fireball matter is affected by the larger k{T) 
values at lower temperatures. 

The total momentum anisotropy e^, on the other hand, 
is more strongly affected by a low-temperature growth of 
k(T) (solid lines in Fig 10b). is suppressed relative to 
the flow anisotropy Sp by the non-equilibrium corrections 
^ TT^'' in the energy-momentum tensor. When the relax- 
ation time T^r is allowed to grow large in the hadronic 
phase, this suppression is found to be reduced, and the 
reduction is relatively larger for large values of (?y/s)HG 
(corresponding to a larger Navier-Stokes value Trfjg) than 
for smaller (?//s)hg- We also note that this suppression 
of the 7r'^''-contribution to is visible already at early 
times when the larger k(T) values affect only the fireball 
corona. In fact, for constant 77/5 = 0.16 (solid squares 
and circles) the low-temperature growth of k(T) leads to 
a bigger increase of over Sp at early than at late times; 
this is due to the larger longitudinal expansion rates at 
early times which lead to larger Navier-Stokes values for 
A = TT^^— TT*'*' everywhere, thus causing greater sensitiv- 
ity to increased k(T) values in the fireball corona. In the 
case of T-dependent 77/s (solid upright and inverted tri- 
angles) the effects from a delayed response are larger at 
late times; in this situation, the Navier-Stokes values for 
A = 7r'^^— TT^^ grow in the hadronic phase due the sudden 
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FIG. 11: (Color online) Differential elliptic flow V2(pt) for 
charged hadrons, using a temperature dependent k(T). Same 
parameters as in Fig. 10. 



increase of 77/5 below Tc, clearly reflected by a "kink" in 
the growth of e'p around t—tq = 4 fm/c (see upright green 
solid triangles in Fig 10b). This kink is largely washed 
out by a simultaneous rise of k{T) (inverted blue solid 
triangles in Fig 10b). 

The behavior of the total momentum anisotropy 
is directly reflected in the charged hadron elliptic flow, 
shown in Fig. 11. We point especially to the reduction 
of the (negative) tt'^^ contributions to in the case of 
T-dcpendent {r]/s){T), which manifests itself through re- 
duced 6f corrections to V2{pt) which again are most pro- 
nounced at large pt (green triangles and blue inverted tri- 
angles in Fig. 11). For constant 77/s, on the other hand, 
the larger hadronic relaxation time has little effect on the 
differential V2{pt), consistent with the very small effect 
on the total momentum anisotropy at late times seen 
in Fig. 10b. 



VII. DISCUSSION AND CONCLUSIONS 

Figure 11 has important implications: Comparing the 
blue line with inverted triangles to the case of constant 
K and r]/s (black squares), we conclude that the suppres- 
sion of V2{pt) reflected in the blue line could have arisen 
in two different ways: (i) by a large increase of i]/ s in the 
hadronic phase, accompanied by a similarly large increase 
of K, as shown here, or (ii) by a much less pronounced 
increase of the hadronic shear viscosity, compensated by 
a correspondingly reduced increase of the hadronic relax- 
ation time. In other words, the hadronic shear viscosi- 
ties and relaxation times extracted from a given charged 
hadron V2(j>t) are strongly correlated and impossible to 
determine independently from a single elliptic flow mea- 
surement. Whether and how the systematic exploration 
of differential elliptic flow for different particle species 
and different collision systems at different centralities can 
help to resolve this ambiguity remains to be seen. 

The study presented here shows that any discussion of 
large dissipative effects in the hadronic phase of heavy- 
ion collisions, reflected by specific shear viscosities and 
(scaled) microscopic relaxation times that grow as the 
system cools below the critical quark-hadron transition 
temperature, is really a discussion of 6f, i.e. of the devia- 
tion of the freeze-out distribution function from its local 
equilibrium form and its reflection in the flnal hadron 
spectra and anisotropics. As the system cools and ap- 
proaches kinetic freeze-out, dissipative effects become 
stronger and stronger, bringing the framework of vis- 
cous hydrodynamics closer and closer to breakdown. In 
this sense, our results have to be taken as qualitative 
insights but should not be confused with quantitative 
predictions. Their main value, as we see it, is that they 
shed light on and help to classify and qualitatively under- 
stand the late-stage dissipative effects on hadron spectra 
and their elliptic flow as seen in a realistic microscopic 
approach (as embodied, for example, by VISHNU). The 
results presented here do provide support to the conclu- 
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sion of Ref. [11] that an effective viscous hydrodynamic 
description of the hadronic stage in heavy-ion colhsions, 
if vahd at ah, hkely requires both large shear viscosity 
and long relaxation times below Tc. 
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